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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL NOTE 2694 

A METHOD FOR STABILIZING SHOCK WAVES IN CHANNEL 
FLOW BY MEANS OF A SURGE CHAMBER^ 


By Stanford E. Neice^ 


SUMMARY 


In order to stabilize normal shock waves in channel flow against 
the effect of distxirbances originating downstream, a method based on 
mass removal from the channel by means of a surge chamber was developed 
and experimentally tested in an intermittent blowdown-type wind tunnel 
at Cornell University. A theoretical analysis of the flow in a channel 
shape similar to that used in a typical double-throat supersonic wind 
tunnel indicated that the mass-removal technique was effective in damping 
the motion of the normal shock caused by a strong compression pulse 
originating downstream. The results of experimental tests were in quan- 
titative agreement with the theoretical analysis. Further experiments 
indicated that the mass-removal technique was effective in damping the 
oscillatory motion of the normal shock caused by continuous small, random, 
downstream distiarbances . 


INTRODUCTION 


An important factor in the design and operation of supersonic wind 
tunnels is their large power consmnption. A large part of this loss 
results from the fact that supersonic channel flow can be converted into 
subsonic flow only through a normal shock wave (see ref. l) . A success- 
ful method for reducing this power loss is to lower the Mach number at 
which the shock occurs. This method has led to the use of the "double- 
throat" type of supersonic wind tunnel. 

For the most efficient operation of the double -throat wind tiuinel, 
the normal shock is placed in a position just downstream of the second 
minimum section. However, distiirbances originating downstream in the 
diffusor and exhaust system will interact with the normal shock and cause 
^The body of this report is a thesis which was submitted in February 1950 
in a partial fulfillment of the requirements for the degree of Master of 
Aeronautical Engineering in the Graduate School of Aeronautical Engineer- 
ing, Cornell University, Ithaca, New York. 

^Now at the Ames Aeronautical Laboratory, Moffett Field, Calif. 
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it to be displaced from its original position. If these disturbances 
are strong enough, they can cause the shock to be permanently displaced 
to a position upstream of the test section. Subsonic flow will then 
exist in the test section and the wind tunnel is said to be "unstarted." 
In order to avoid this highly undesirable condition, it is often neces- 
sary to place the normal shock quite far downstream of the second mini- 
mum section. Operation with the shock far downstream of the second 
throat results in an Increased power loss and defeats partially the pur- 
pose of the second minimum section. 

The power consumption of the double -throat wind tunnel could be 
reduced if there were a suitable method for stabilizing the normal shock 
against the effect of these disturbances, thus permitting the shock to 
be placed closer to the second minimum section. It is the purpose of 
this paper to show how the use of a surge chamber, connected to the chan- 
nel immediately upstream of the operating position of the normal shock, 
can produce the desired stability. 

The author is grateful to Professor Arthur R. Kantrowitz for his 
suggestion of the topic, advice throughout the investigation, and for the 
information contained in Appendixes A and B. 


SYMBOLS 


area of transverse slit 
speed of sound 

specific heat at constant pressure 
specific heat at constant volume 
Mach number 
mass flow 

characteristic quantity, 
pressure 

characteristic quantity, 
gas constant 
area of channel 
temperature 
time 
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U velocity of shock wave 

u velocity of fluid 

X station along channel 

7 ratio of specific heats 

P density 

T entropy 

Subscripts 

a conditions behind steady-flow normal shock 

b conditions behind leading edge of pulse 

c conditions at rear of pulse 

s conditions relative to moving shock wave 

o standard conditions 

1,2, 3^ 4^ 5 conditions in various sectors of the characteristic 
diagram 


Superscript 


average values 

THEORETICAL CONSIDERATIONS 

Stability of Normal Shock Waves 
in Channel Flow 


In the normal operating condition of a conventional double -throat 
supersonic wind tunnel, the shock wave converting the supersonic flow 
to subsonic flow will occur downstream of the second minimum section as 
shown in sketch (a). The normal shock wave in this case is in its pri- 
mary stable equilibrium position which means that it will return to this 
position after undergoing small displacements. 
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Sketch (a) 

There are two other equilibrium positions, also shown in sketch (a), 
which occur at points in the channel where the channel area is the same 
as the area at the primary stable equilibriijm position. At the indicated 
position in the converging part of the channel, the normal shock is in a 
condition of unstable equilibrivun. Consequently, if the shock were to 
exist in this position, then any displacement, however small, would cause 
the shock to assume either stable position, depending upon the direction 
of the displacement. 

Displacement of the normal shock from the primary equilibrium posi- 
tion can be caused by disturbances originating downstream. These dis- 
turbances propagate upstream through the subsonic flow in the diffusor 
section and interact with the main shock. The resultant shock is then 
set in motion and will undergo a displacement proportional to the 
strength of the disturbance. As pointed out in reference 1 , a compres- 
sion disturbeince will cause the shock to be displaced upstream in the 
direction of the unstable equilibrium position. If the compression dis- 
turbance is strong enough to cause a displacement of the shock upstream 
of the unstable equilibriimi position, then, as has been pointed out, the 
shock will assume the secondary stable equilibriiam position shown in 
sketch (a). Subsonic flow will then exist in the test section, and the 
tunnel will be unstarted. 


Operation of the Surge Chamber 


In the method presented herein, the shock position is stabilized by 
the appropriate use of a surge chamber, which, in effect, creates a 
coianterdisturbance which neutralizes the action of a compression disturb 
ance coming from downstream. 


MCA TW 2694 


5 


The surge chamber is connected to the channel (see sketch (b)) by 
transverse slits at a point just upstream of the primary equilibrium 
position. In the normal running condition, the surge chamber is at the 
same pressure as that existing in the channel at the position of the 
slits. The normal shock, set in motion by interaction with a compression 
disturbance, crosses the slits, thus raising the pressure in the channel, 
and causes air to flow out of the channel into the surge chamber. 



Sketch (b) 

The removal of mass in this manner, in effect, produces an expansion 
which propagates throughout the channel. Part of this expansion inter- 
acts with the upstream-moving shock wave and dissipates some of the 
energy supplied by the initial compression distxirbance . By varying the 
opening of the transverse slits, a sufficiently high rate of mass 
removal may be obtained to bring the shock wave to rest before it passes 
the unstable equilibriiim position. A normal shock brought to rest in 
this manner will then return to its primary stable equilibrium position 
downstream of the second mi m' mum section. 

By the application of the method of characteristics for nonsteady 
one-dimensional flow, it is possible to trace the interaction of the 
normal shock and compression pulse on a time -displacement diagram and to 
determine the effect of mass removal to the surge chamber. In order to 
analyze this problem, however, it was necessary to develop an extension 
to the existing method of characteristics which includes the effects of 
mass removal from a point in the channel. 
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The Method of Cheiracterlstlcs Including the 
Effects of Mass Removal at a Point in the Channel 


The method of characteristics for isentropic, inviscid, nonsteady, 
one -dimensional flow in a channel of constant cross section was first 
introduced in l859 by Riemann (ref. 2). In reference 3, the work of 
Riemann was extended to include variations in entropy and channel area, 
while in reference 4 the effects of heat transfer were considered in 
some detail. The general procedure followed in those papers, and 
exploited in the present analysis, was to define certain characteristic 
quantities from considerations of the one-dimensional equations of 
motion and continuity. These characteristic quantities, designated as 
P and Q,® are associated with the two families of characteristics in the 
flow field. 

In the present analysis the equations of continuity and motion are 
modified to include mass removal from the channel, and the effect of 
mass removal on the characteristic quantities, P and Q, is determined. 

The equation of continuity for one -dimensional flow is 

S ^ ^ (puS) = 0 ( 1 ) 


If, however, mass is removed from the channel through openings in 
the walls, the continuity equation takes the following form: 


„ dp d /„ dm „ 


( 2 ) 


In the present analysis it is assumed that the mass is removed from 
the channel at a constant rate ( ^ = 0 ) . Consequently, the equation 

\dt J 

of motion for one -dimensional flow remains unchanged as 


^+u^+i^=0 (3) 

8t dx P dx 


^These quantities, as used in references 3 and 4 as well as the present 
analysis, differ only by a factor of 2 from the Riemann invariants as 
given in reference 
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With the aid of the first law of thermodynamics, the equation for 
the entropy of a perfect gas can be written in the following differential 

T DTTn R • 


~ = — d(loge a) - d(loge p) 
dT) 2cp 

^ = — <i(loge a) - d(loge p) 




By the use of equations ( 4 ) and 
gas in the form p/ p = rt^ equations 
following forms ^ respectively ; 


the equation of state for a perfect 
(2) and (3) can be brought into the 


2cy Sa 2c V Su ^(log^ S) 

— + u — + a — + ua 1 

R St R bx Sx 


a 

R ^^t 


+ 



a 

pS hx 


0 


\ 


> 

J 


( 5 ) 


Bu 2c V ^a a^ 

— + u — + — 1 a 0 

ot dx R dx 7R dx 


( 6 ) 


By adding and subtracting equation (6) from equation (5), the 
following relations are obtained: 


d 

dt 


d 

3 t 




+ (u + a) — ( — 4 a + u + ua 

V R / ax 

7R dx 
+ (u 

aD 4 

7 R dx " R Dt PS dx " ° 


a ^ ^ ^ _ 

R Dt ^ d^ " ° 


r) ^ „ A ^(log„ S) 

a; I — 1 a - u I + ua ^ 

dx V R / ^x 


\ 



J 
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where the notation 


= i. + u -i 

Dt 8t dx 


( 8 ) 


defines an observed rate of change while moving with the fluid particles 
at a velocity u. 

The characteristic q.uantities, P and Q, are defined as 


2cy 

P = a + u 

R 

2Cy 

Q = — - a - u 
R 


( 9 ) 


and the observed rate of change while moving along a characteristic line 
in the fluid is defined symbolically as 


— = ^ + (u ±a) ( 10 ) 

6 t 8 t ox 


where the plus (+) or minus (-) signs are associated with P and Q fam- 
ilies of characteristics, respectively. 


From the definitions of equations ( 9 ) ajid (lO), equation ( 7 ) can 
be written as 


6P _ 8(logeS) ^ ^ a DT) _ a ^ 

5t 8x yR 8x R Dt p3 dx 


= -ua 8 (loge S) + 

6t Sx yR dx R Dt pS dx 


( 11 ) 


which expresses the rate of change of the characteristic quantities 
along their respective paths. 
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From equation (lO), the following can be deduced: 


Along a P characteristic 



dx 


Along a Q characteristic 


( 12 ) 


6t 

Q 



9 


The channel area of a wind tunnel will depend only on the station 
along the channel. Consequently, the rate of change of this quantity 
along P and Q characteristics, with the aid of equation (lO), is 


Since the removal of mass from the channel has been assumed to take 
place at a constant rate, the change in mass flow is, therefore, a 
fimction only of the distance along the channel. From equation (lO) the 
rate of change of this quantity along the characteristic lines can be 
expressed as 


S(loge S) 
6t 


P 



(13) 


5 (loge S ) 
&t 


= (u - a) 


Q 


5m 

6t 


= (u + a) ^ 


P 


( 1 ^) 


5m 

5t 



Q 


j 
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With the aid of equations (12)^ (13) ^ and (l^), and noting that the 
particle entropy change Dil/Dt will be zero In the absence of heat 
transfer and shock waves, the variations In the characteristic quantities 
will then be 


ua 


5P 5(loge S) + — 6T1 — 6m 


u + a 


7R 


pS(u + a) 


6Q = - 


u a 


u - a 


6(loge S) + — 6T1 

7R 


pS(u - a) 


6m 


(15) 


It should be pointed out that the quantity 6m will Itself have a 
positive or negative sign, depending upon the direction (upstream or 
downstream) which a characteristic crosses the point of mass removal. 
The quantity 6m can also be expressed in the following form: 


C ^ A 

6m = — Ax 
ox 


(16) 


where 6m/6x is positive if mass is removed from the channel and Ax 
will be positive for the downstream direction. 


The magnitude of the quantity 6m can readily be evaluated by con- 
sideration of the equations of continuity and momentum referred to the 
dimension perpendicular to the channel. A simplification is possible, 
however, since the pressure rise across the shock wave, the pressure 
difference between the channel and the surge chamber, is generally 
sufficiently large so that the flow through the slits is choked. In 
this case, therefore, 

1 6m I = p* a* A (1?) 

and the variation of P and Q, produced by mass removal only, reduces to 


I6P[ = 


p a A 


(u + a) P a S 


5q! = 


* A 

P a A 


(u - a) p a S 


(18) 
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In equation (18), the ambient stream quantities in the channel correspond 
to stagnation conditions with respect to the cross flow through the 
transverse slits. Consequently, the following ratios are constants: 


P* 

P 


0.634 


(19) 


^ = 0.913 

cL 


The resultant magnitudes of P and Q, due to mass removal, and assigning 
the flow through the slits is choked, are 


|8P| = 0.578 


a^ A 
(u + a) S 


l5Q] = 0.578 - 

(u - a) S 




( 20 ) 


The effects of shock waves and large entropy discontinuities are 
not included in the foregoing theoretical considerations. Both these 
effects were present in the analysis by the method of characteristics. 

A discussion of how these effects are handled is given in Appendix A. 

In performing the characteristics analysis, it was also necessary 
to consider three types of interactions involving shock waves; namely, 

1. Interaction of the pulse shock with the steady-flow normal 

shock 

2. Interaction of the resultant shock with the flow out of the 

channel at the instant the shock crosses the transverse 
slits which open into the surge chamber 

3. Interaction of the resultant shock with expansion waves 

from downstream which "catch up" with the shock and 
change its strength 

The method of solution for each of these interactions consists of a 
series of successive approximations and is discussed in Appendix B. 
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APPLICATION TO PARTICULAR DESIGN 


Preliminary Considerations 


Shape of channel s- The method of characteristics as developed herein 
was used to study the stability of the steady-flow normal shock in a 
small double-throat supersonic wind tunnel equipped with a surge chamber. 
The contour of the channel is given in figure 1. The Mach numbers at the 
test section and the second throat were 1.64 and 1.42, respectively. The 
supersonic part of the channel was designed according to the steady-flow, 
two-dimensional method of characteristics. A diffusor half angle of 2P 
was used in order to minimize separation losses in the subsonic part of 
the flow. In the operating condition, the equilibrium position of the 
normal shock occurred 5*8 inches downstream from the first minimum 
section, that is, I .3 inches downstream of the second throat and O .3 
inch downstream of the transverse slits. The Mach number immediately 
upstream of the normal shock was 

Construction of velocity profiles .- With the shape of the channel 
and the position of the normal shock now specified, the steady-flow 
velocity profile was computed and is shown in figure 2. The initial 
properties of a pulse originating downstream in the diffusor section are 
generally not known since they may arise from a variety of causes. It 
was shown in reference 1, however, that the velocity profile of a pulse 
takes a definite shape as it travels upstream, regardless of the condi- 
tions at its origin; the upstream section of a compression pulse will 
steepen to form a leading-edge shock, while the velocity profile of the 
expansion phase of the pulse behind this shock will assume a slope given 

by 


All the disturbance pulses considered in this analysis are compres- 
sion pulses having this shape of velocity profile. Two such pulses are 
shown in figure 2 , superimposed on the steady-flow velocity profile at 
the position of impending intersection with the steady-flow normal shock. 

Determination of initial flow conditions .- In order to start the 
characteristics analysis, the flow conditions, including the values of 
the characteristic quantities (eq. ( 9 )) must be determined at all points 
in the channel at the time of the intersection of the pulse shock with 
the steady-flow normal shock. The details of this computation are given 
in Appendix C . 



pulse 



flow 
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Analysis and Results 


The method of characteristics was first applied to the channel, 
without mass removal, to find the "critical" pulse, that pulse which will 
supply sufficient energy to displace the normal shock to the unstable 
equilibrivun position and, hence, unstart the flow. The order of magni- 
tude of the critical pulse can be determined according to the methods 
of reference 1 where it was shown that, in the vicinity of Mach number 1, 
the velocity-profile area of the critical pulse'*' will be equal to the 
area between the steady-flow velocity profiles formed by the normal 
shock in the stable and unstable equilibrium positions (see fig. 2). 

This rule will, of course, give errors for Mach numbers greater than 1 
but can be used as a first approximation. The present analysis, using 
the method of characteristics, showed that the critical pulse had a 
velocity-profile area about 68 percent of that given by reference 1. 

This result is in agreement with an investigation described in reference 
6 which indicated a velocity -profile area of approximately 65 percent. 

The critical pulse is the smaller pulse shown in figure 2. A time- 
displacement history of the resultant shock formed by the interaction of 
the critical pulse with the steady-flow normal shock is given in figure 3. 

It was decided, as an appropriate measure of the effectiveness of 
the surge chamber, to determine the size of transverse slit necessary 
to prevent unstarting of the flow by a pulse having a velocity-profile 
area double that of the critical pvilse already determined. This double 
critical pulse is the larger pulse shown in figure 2. Figure 3 also 
presents time -displacement diagrams for the resultant shocks formed by 
the Interaction of the steady-flow normal shock with the double critical 
pulse for the following conditions: no mass removal, l/8-inch transverse 

slits, l/U-inch transverse slits. It is evident that l/4- inch slits are 
close to the smallest size which will dissipate the action of the double 
critical pulse. 

In all cases, the pressure in the channel, at the point of mass 
removal, was sufficiently large to cause the mass flow into the surge 
chamber to be choked, thereby allowing the use of equation (19) to com- 
pute 5P and 5Q across the slits. 


EXPERIMENT 


Part of the experimental phase of this investigation consists in 
generating strong compression pulses in the flow and observing the inter- 

action of the pulses with the steady-flow normal shock for the wind 

^The velocity-profile area of a pulse is the area between the steady-flow 
velocity profile of the channel and the velocity profile imposed by the 
pulse . 
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tunnel in one of the two following operating conditions: (a) slits into 

surge chamber closed; no mass removal; and (b) slits into surge chamber 
open l/k inch. By means of a modified shock tube, it was possible to 
generate pulses having a velocity profile approximately the same as 
those used in the characteristics analysis. Hence; these experimental 
tests could parallel the analytical studies by determining the relative 
strengths of the strong pulses which would just unstart the flow under 
the stated operating conditions. 

A further check on the effects of mass removal was made by operating 
the wind tunnel with the normal shock just downstream of the transverse 
slits and observing the effectiveness of the surge chamber in damping 
the movement of the shock caused by the incidence of small random dis- 
turbances. Since small random disturbances are encountered more often 
than strong pulses; it was felt that this test might be more representa- 
tive of actual operating conditions. 


Apparatus 


Wind tunnel and optical system .- The experimental tests were car- 
ried out in a small intermittent supersonic wind tunnel in the gas dynam- 
ics laboratory of the Graduate School of Aeronautical Engineering; 

Cornell University. The width of the channel was 1 inch. The only modi- 
fications to the tunnel were the construction of new nozzle blocks; as 
shown in figure 1 ; and the construction of a bracket for mounting the 
pulse -generating device at the end of the diffusor. The transverse 
slitS; which connect the channel to the surge chamber; were adjustable 
to allow an opening up to 1/2 inch. The surge chamber was designed to 
be large enough so that the pressure would not change appreciably within 
the length of time necessary for the shock to travel from the stable 
equilibrium position to the unstable equilibrium position and return. 

The optical apparatus consisted of a schlieren system of conven- 
tional design and was provided with two light sources. A continuous 
light source was used for visual observation and alsO; in conjunction 
with a conventional camera shutter; for taking schlieren pictures with 
exposures of 1/25 second. The other light source; used for photography 
only; was a high-intensity spark having an effective duration of 1 micro- 
second. 

Pulse tube.- The compression pulses were generated in a small shock 
tube which was moiinted on the end of the diffusor. The shock tube; 
called a "pulse tube;" differs from the conventional design; described 
in reference 7 ; in that the rearward wall of the high-pressure chamber 
is placed close to the diaphragm. The longitudinal dimensions are pro- 
portioned so that the front of the expansion wave reflected from the back 
wall just overtakes the shock at the forward end of the tube. At this 
instant the velocity distribution over the length of the tube decreases 
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aljnost linearly from a large negative value behind the shock to a very 
small positive value at the rear of the pulse. Figure 4 presents a 
typical analysis by the method of characteristics of the flow in the 
pulse tube together with a plot of the velocity profile. It will be 
recalled that the velocity profile of the compression pulses used in the 
theoretical analysis (see fig. 2) jumps discontinuously through a shock 
wave to a large negative value and then increases almost linearly to 
zero. Hence, the pulse tube generates a compression pulse with the cor- 
rect shape of velocity profile needed to check the prediction of the 
theory. 

The following dimensions of the pulse velocity profile are required: 
(a) the velocity discontinuity at the shock front, and (b) the velocity- 
profile area. The velocity discontinuity at the shock front of the pulse 
is determined by the initial pressure ratio across the diaphragm. The 
area of the pulse can be adjusted to the required value by changing the 
over-all length of the pulse tube. The pulse tube can be designed to 
give precisely the required velocity-profile area corresponding to only 
one particular value of the velocity discontinuity at the shock front. 

At other values of velocity discontinuity, the profile area obtained in 
the pulse tube will differ slightly from the profile area required in 
the theoretical analysis. This difference, however, was found to be less 
than 4 percent throughout the range of values employed in the experiment. 

The pulse tube was designed according to the dimensions of the 
double critical pulse shown in figure 2. This pulse had a velocity dis- 
continuity, Au/uq, across the shock front of O.515 and a velocity -profile 
area of 1.57 (units of u/ao X inches). It was anticipated, however, 
that there would be losses in pulse strength caused by the pulse travers- 
ing the gap between the end of the pulse tube and the end of the dif- 
fusor, as well as losses resulting from the separated flow in the 
diffusor. It was estimated that these losses would cause a reduction in 
the over-all pulse velocity profile of 20 percent. For piirposes of 
design, therefore, the dimensions of the double critical pulse were 
increased by 20 percent which brought the velocity-profile area to a 
value of 1.90 (units of u/aQ x inches) and a leading-edge velocity 
discontinuity, Au/uq, of O.6IO. Two additional analyses were carried 
out for weaker pulses. From these analyses the relation was determined 
between the velocity-profile area and initial pressure in the high- 
pressure chamber of the pulse tube. The graph of this relation is pre- 
sented in figure 5. 

Figure 6 shows the entire test setup with the pulse tube in its 
operating position at the end of the channel. A larger photograph of the 
pulse tube is shown in figure 7 • Also shown in figure 7 is the diaphragm 
rupture mechanism which consisted of a needle attached to a solenoid 
switch. Closing the circuit of the solenoid switch caused the needle to 
be displaced inward, thus rupturing the diaphragm. Due to the deflection 
of the diaphragm under pressure, a filler plate was attached to the back 
wall of the high-pressure chamber (see fig. 8). The size of this plate 
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was approximately equal to the added volxme caused by the diaphragm 
deflection and^ consequently, acted to keep the pulse size closer to 
that dictated by the theoretical design. 


Procedure 


In order to study the effect of mass removal against the action of 
the compression pulses from the pulse tube, tests were conducted with 
the transverse slits l/h inch open and were repeated with slits closed. 
The procedure in both cases was to fill the high-pressure chamber in the 
pulse tube to a desired pressirre, then start the wind tunnel and bring 
the normal shock to its proper position downstream of the second throat. 
The diaphragm of the pvilse tube was then pxmctured. The effect of the 
interaction of the compression pulse with the steady-flow normal shock 
was observed on the schlieren screen: Either the normal shock returned 

to its stable equilibrium position downstream of the second throat or was 
driven to the secondary stable equilibrium position upstream of the test 
section. These tests were repeated with different values of pulse-tube 
pressure \mtil that pressure corresponding to the critical pulse was 
determined. 

In order to determine the effect of mass removal against the action 
of RTtiai 1 random disturbances, a second series of tests was carried out 
with the air flow from the end of the diffusor partially blocked. The 
sharp diversion of air at this point created random pressure variations 
which were propagated upstream through the diffusor. It was felt that 
the disturbance level created in this manner would correspond roughly 
to disturbances produced by the exhaust system of an actual wind tunnel. 
For this part of the experiment, the normal shock was placed in a posi- 
tion at the downstream edge of the transverse slits. The effect of mass 
removal to the surge chamber was observed by comparing the oscillation 
of the normal shock when the slits were l/k inch open with the oscil- 
lations present when the slits were closed. An attempt was made to show 
this effect by taking schlieren photographs with a spark ( 1 -microsecond 
exposure) and with an exposure of I/25 second. The normal shock wave 
was photographed at both these exposures, first with the slits open 
l/h inch and then with the slits closed. If mass removal is effective 
in damping oscillations of the shock wave, there will be a close simi- 
larity between the l/25-second and 1-microsecond exposures. 


RESULTS AND DISCUSSION 


The experimental tests to determine the effect of mass removal 
against the action of strong compression pulses revealed that the pulse 
necessary to unstart the tunnel with l/ 4 -inch transverse slits was 
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produced by a pulse-tube-chamber pressure of approximately 62 psi gage. 
With the transverse slits closed, the pulse necessary to unstart the 
tunnel was produced by a pressure of about 28 psi gage. As shown in 
figure 5^ the velocity-profile areas corresponding to these pressures 
are I.80 and O.87, respectively. The velocity -profile area of the pulse 
required to unstart the tunnel with l/ 4 -inch slits is approximately 
double the area required with slits closed. This result is in agreement 
with the theoretical analysis, as shown in figure 3, which predicts that 
the action of the double critical pulse will be dissipated by mass 
removal through l/ 4 -inch transverse slits. 

The second series of tests demonstrated the effect of mass removal 
against the action of small random disturbances. With the slits open 
1/4 inch, the normal shock was never driven upstream by the random dis- 
tuxbances, and a marked decrease in the oscillation of the normal shock 
was noted as compared with the oscillations present when the slits were 
closed. 

ScMieren photographs, with exposure times of I/25 second and 1 
microsecond, were taken with the transverse slits both closed and open 
and are presented in figures 9 and 10 , respectively. The oscillatory 
motion of the normal shock, caused by the random disturbances, can be 
observed in figure 9 by noting the difference in clarity between the 
photographs taken at the two exposure times. The effect of mass removal 
in damping these oscillations is to be observed by comparing the relative 
clarity between the l/ 25 -second exposures of figures 9 and 10 and, also, 
by noting the similarity between the l/25-second and 1-microsecond expo- 
sures of figure 10. 

When taking schlieren photographs, the shock could not be positioned 
visually because the visual screen was blocked by the photographic plate 
holder. An approximate shock position was fixed by reference to the 
settling chamber pressure. Under these conditions it was difficult to 
photograph the normal shock at a desired average position. As a result 
of this difficulty, the photographs of figure 10 show the normal shock 
at an average position farther downstream from the transverse slits than 
was desired. Although some damping is Indicated by the relative clarity 
of the l/25-second exposure in figure 10, as compared to the 1-microsecond 
exposure, the effect is not as striking as was observed visually. 


CONCLUDING REMARKS 


In order to stabilize normal shock waves in channel flow against 
the effect of disturbances originating downstream, a method based on 
mass removal from the channel by means of a surge chamber was developed 
at Cornell Unlvers'tv. The stabilizing action is initiated by the motion 

of the ■':'='ac'y* ' w normau shov.k caused by the disturbances and requires 
the u. no add.' clonal power facilities. 
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With the use of the method of characteristics (extended to include 
mass removal at a point in the channel), the flow was analyzed in a chan- 
nel shape similar to that used in a double-throat supersonic wind tunnel. 
It was found that the mass -removal technique was effective in damping the 
motion of the normal shock wave produced by strong compression disturb- 
ances . 

To check the findings of the theoretical analysis^ a series of 
experimental tests was performed in a small double -throat wind tunnel into 
which compression disturbances of known strength were introduced by a 
pulse generating device. The experimental results were in quantitative 
agreement with the results of the theoretical analysis. An additional 
series of experimental tests was performed in which it was fo\md that the 
mass -removal technique was effective in damping the oscillatory motion 
of the steady-flow normal shock caused by small random disturbances. 

From the theoretical and experimental results, it can be concluded, 
therefore, that the mass -removal technique is effective in stabilizing 
normal shock waves against compression disturbances originating down- 
stream, and that the performance can be accurately computed by the method 
of characteristics. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 31^ 1953 
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APPENDIX A 

APPLICATION OF THE METHOD OF CHARACTERISTICS 
ACROSS STRONG DISCONTINUITIES 
Normal Shock Waves 


The change in each of the characteristic quantities, P and Q, and 
the entropy, t| , across a normal shock wave may be computed using the 
customary Rankine-Hugoniot relations. However, in the construction of 
characteristic diagrams where solutions across shock waves are required 
many times, it has proved convenient to prepare two auxiliary charts to 
facilitate these computations. To begin with, it is desirable to define 
the two classes of shock waves. A P shock is defined simply as a shock 
wave which, on a time-displacement diagram, is moving to the right rel- 
ative to the fluid. Conversely, a Q shock is a shock wave that moves 
to the left relative to the fluid. In addition, all velocities are taken 
to be positive for movement from left to right and negative from right 
to left. 


The first of these auxiliary 
charts expresses the change across 
a shock wave in the quantity P in 
terms of the corresponding change 
in Q. Consider a P shock in 
channel flow as shown in sketch (Al) . 
The changes in the characteristic 
quantities, expressed in non- 
dimensional form by dividing by the 
speed of sound, ai, is 


2 

-^u ' 



U 2 



Sketch (Al) 



2c 


2^ 

R 


ai - ui 


ai 


'^Qq - Qi^ 


2Cy 2Cy 

^ — ag - U2 - ~ ai + ui 


ai 


> (Al) 


J 
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-(U-Ug) 

-(U-u,) 
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Upon transferring to a coordinate 
system moving with the normal shocks 
as shown in sketch (A 2 ), it is found 
that the relations given in equation 
(A1) now become 


Sketch (A2) 

2.Cy 


■ Pi 

ai ✓ p 


R 


82 - (u - Ug) - ^ ai + (U - Ui) 

ai 


> (A2) 




2 c, 


2 c, 


^ 83 + (U - U 2 ) - ~ ai - (U - ui) 


ai 


y 


Equation (A2) is essentially unchanged from equation (Al) since it merely 
involves the addition and subtraction of the shock velocity, U, in the 
numerator of the right side. This observation serves to illustrate the 
fact that (P 2 -Pi)/ai and (Oe-Qi)/ai are independent of coordinate 
system. Equation (A2) may be brought into the following form; 



/ (A3) 


( 


Qg 




In terms of quantities listed in the standard steady-flow normal - 


shock table, 


M 


M 


U - Ui 

“ ai 


m 
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Equations (A3) then becomes 



(A5) 



With reference to any of the standard steady-flow normal-shock tables, 
the chart of (Pa-Pi)/^! as a function of (Oe-Qi)/ai can be constructed. 


For comparison, a Q shock in 
channel flow would be represented as 
in sketch (A3) . Upon transformation 
to a coordinate system moving with 
the shock velocity, U, the flow would 
be represented as in sketch (a4) . 

The change in the characteristic 
quantities for the case of the Q 
shock would now be 


/ 

U -1 

2 


»2 


Sketch (A3) 


/ 

2 

ITTTTTTTTTTTTi 

(-u*u^r 

\nnmiinr 


Sketch (a4) 



^2 /2c^ -U + U2\ /2c^ -U + ui\ 

ai V~ ag ) ■ ) 


(A6) 


ai 4 


ag f 2c^ -U + U2^ 

ai V R J 



-U + ux' 


0,1 
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In terras of quantities listed in standard steady-flow normal-shock 
tables 


M 


Si 


-U + ui 
ax 



-U + Ug 
^2 






and equations (A6) then become 



(A7) 


(A8) 


Comparing equations (A8) with the corresponding equation for the 
P shock, the following rule can be stated: 


( 


Pa 



Qq ~ Qi 
ai 


Q 


Qg - 
ai 


P2 - P f 

ai 


Q 






y 


(A9) 


The curves of (P 2 “Pi)/^i a function of ^'or a P 

shock can^ therefore^ be used for Q shocks by merely interchanging the 
ordinate and abscissa. 
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The second of the necessary charts concerns the entropy change 
across a normal shock. The following relation expresses the entropy 
change between any two eq.uilibrii;m conditions of a perfect gas. 



7R 7R 




(AlO) 


For a normal shock, however, the ratios a^/ai and Pg/pi are related 
through the Rankine-Hugoniot equations. With the use of a standard set 
of steady-flow normal-shock tables, it is therefore possible to compute 
(’l2“^l)/7R as a fimction of only. 


In applying the method of char- 
acteristics across a normal shock 
wave, a complete solution involves a 
knowledge of the six quantities as 
shown in sketch (A5) . In the present 
analysis, which is concerned only 
with shocks, the quantities Pi, 

known, either 

from boundary conditions or from 
previous computations. With the aid 
of the chart of ( P 2 -Pi ) /ai as a 
function of (Q 2 "^i)Ai> "the value of 
Pg can be determined. According to 
the definitions of P and Q, the flow 
velocity and speed of sound may be 
computed from the following general 
relation: 



Sketch (A5) 


P + Q 
^c^/r 


> 

P - Q 

i 


(All) 


u 


2 
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With the value of a^/ai now known, can be determined from the chart 

of (T)2“’ni)/7R as a function of Sg/ai* 

Strong Entropy Discontinuities 

When a large entropy discontinuity exists in a fluids the pressures 
and velocities on each side of the interface are identical. For this 
condition, equation (AlO) becomes 


^2 - ^1 n 

— = logp 

7R 7R 



(A 12 ) 


from which the speed-of-soimd ratio between conditions on each side of 
the interface is 


ai 



(AI3) 


Since the velocities on either side of the interface are constant, the 
difference between the characteristic quantity P across the entropy 
discontinuity is expressed as 


2cv , > 2cv 

Pi = — (as - ai) = — ai 


_ 

ai 


(A 14 ) 


Substitution of the expression for the speed-of-soiond ratio in equation 
(AI3) into equation (Al 4 ) gives the following relation: 


P2 - Pi 


2 c, 


ai 


^ 2 - ^1 

2 Cp 


R 


e 


-1, 


(AI5) 
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By a procediire analogous to that used in obtaining equation (AI5), 
the following expression can be written for conditions across the inter- 
face: 


Pi + Qe - 2 — ai 



(AI6) 


Combining equations (AI 5 ) and (A16) yields 


P2 ” Pi - (Pi + Qs) 


% - \ 


e _ 1 


- ^1 
. 2Cp 


+ 1 


or 


T] _ T] 

Pg - Pi = (Pi + Q2) tanh 

4 c r, 


(AIT) 


By a similar procedure to that used in obtaining equation (AI7), the 
following expression can be obtained: 


Qs - Qi = (Qi + P2) tanh (A18) 

4Cp 

For very small entropy discontinuities, the first term of the series 
expansion for the hyperbolic tangent may be sufficient, and equations 
(AI7) and (A 18 ) reduce to the form given in equation (15) in the text of 
this report. 


The application of equations g f 

(AI7) and (A 18 ) to the determination ^ 
of the flow quantities across a typi- 
cal strong entropy discontinuity can 
be shown with reference to sketch (A 6 ) . 
The quantities Pi, Qg, T)i, and are 
known from boundary conditions . The 
problem is to determine Qi and Pg . 

From equations (AI7) and (AI8), the 
values of (Pg - Pi) and (^^g - ^1) can 
be computed from known conditions. Qi 
and P2 can then be computed since 
Q2 Pi 3 -re known. 



Sketch (A6) 
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APPEKDIX B 
TYPICAL INTERACTIONS 


The theory of wave interactions is discussed in reference 5 - How- 
ever^ the solution of any actual interaction process requires a procedure 
for numerical computation, and this aspect of the problem is hardly 
touched upon elsewhere in the literature. Accordingly, an attempt is 
made here to present a practical method of solution for those interac- 
tions encountered in the present analyses. 


Interaction of Two Q Shock Waves 


This interaction is encountered when the shock wave at the leading 
edge of the compression pulse encounters the steady-flow normal shock. 

A schematic diagram of this interaction is given in sketch (B 1 ) . The 

. . . . interaction of the two 

discontinuity 




/Entropy 


Com- 

bined 

shock 


Steady 

ftP^ y 
shock - 



Reflected 

disturbance 


Pulse shack 


X 

Sketch (El) 


shock waves results in 

(1) a single strong 
shock wave which moves 
upstream with a velocity 
intermediate to the 
interacting shocks; 

(2) a strong entropy 
discontinuity which 
moves downstream with 
the velocity of the 
fluid; and (3) a weak 
reflected disturbance 


which, as shora in reference 5^ v-'iH be a rarefaction wave, 
for solution is given in the following outline form: 


The method 


Given: 


1 . All conditions (P, Q, u, a, t) ) in sectors 1 , k, 5. 

2 . Velocity is constant across entropy discontinuity: u^ = u^ . 

3. Q and T are constant across the reflected disturbance. 

Q3 = 0 , 4 } fj 3 = ^ 4 ' 

Find : 

1. All conditions (P, Q, u, a, ti ) in sectors 2 and 3. 

2. The velocity of the combined shock. 
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Auxiliary Charts (See Appendix A) 

!• (Pa“Pi)/a-i vs (^-Qi)/ai for Q shock. 

2. (-^ 2 -\)/ 7 R vs — for normal shock. 

3-1 

3. Steady-flow normal-shock tables. 

Procedure : 


1. Assume (Q2-Qi)/ai. 

2. Obtain (P2-Pi)/ai from charts of (P2-Pi)/ai vs (Q2-Qi)/ai for 

Q shock. 


3. Compute P2, Q2, U.2, 8.3. 

14-. Obtain ^3 from chart of (i12-t]^)/7R vs — for normal shock. 

3-1 

5. Compute (Qs - Q2) from equation (A18) . In equation ^Al8) replace 

subscript 1 by subscript 3 , and use the given value of Q3 
in computing the factor ( Qg + Pg ) . 

6. Compare Qg (computed) with Qg (given) and repeat steps 1 through 

5 until these quantities are equal. 

7. Compute Pg from equation (AI7) replacing the subscripts as 

indicated in step 5* 

8. With Qg from step 6 and P3 from step 'J, compute Ug and 3g 

with the use of equation (All). 


With conditions in sectors 1 and 
2 now known, the velocity of the result- 
ant shock can now be computed. Consider 
the motion of this shock as shown in 
sketch (B2) . When transformed to a 
coordinate system moving with the veloc- 
ity of the shock, the picture then 
assiimes the form shown in sketch (B3) . 
The Rankine-Hugonlot relations may be 
used with reference to this latter 
sketch. As such. 


U = ui - ai 

U = ^2 - 


(Bl) 


Z////////////////////// ////// 


/ 

a 

2 


^2 

/ / U /// U ! ! ! n ! ! / ! ! ! H ! ! n n 

Sketch (B2) 

////////////// / // // // ////// 



Uf-U 

U2~U 
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Sketch (33) 
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The velocity of the shock U may be determined from either of the eq_ua- 
tions (B1) with the aid of a set of steady-flow normal-shock tables. 


A manerical example of this procedure may help. In practice it is 
found convenient to express all velocities, and also P and Q,, as multi- 
ples of a standard velocity. In this case the standard velocity is taken 
as the standard speed of soxmd, aQ (1120 ft/sec), and the values of u, 
a, P, and Q listed are actually u/aQ, a/a^, pAq, and qAq. The follow- 
ing initial conditions are given in sectors 1, 4, and 5‘- 


Pi = 5*384 

Ul 

= 1.274 

Qi = 2.836 

ai 

= 0.822 

Ps = 5*384 

^5 

= 0.654 

Q5 = 4.131 


= 0.957 

P 4 = 5*450 

^4 

= 0.290 

Q4 = 4.870 

^4 

= 1.034 


7R 


= 0 


% 

7R 


0.069 


7R 


0.079 


7R 7R 

A ~ ^4 ~ 4.870 

Following the specified procedure, the resulting solution is 

1. Assume (Q2-Qi)/ai = 2.648 

2. Obtain (P2-Pi)/ai - 0.227 from chart of (Q2-Qi)/ai vs (P2-Pi)/ai 

3. Compute P2 = 5*571 U2 = 0.280 

Q2 = 5.011 a^ = 1.058 

4. ag/ai = 1.288, hence (Tl2"'ni)/7R /rR = 0.212 

from chart of Ar,/7R vs a^/ai 
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5 * ( VVAr = 0*133 (n3-Tis)/^p = 0.0133 

tanh (Ti3-T)2)/4cp = 0.0133 

Q3 + Pg = 10.441 Use Q3 (given) 

Q2 - Q3 = 0.139 

6. Qs = 4.871 (computed) as compared to 
Q3 = 4.870 (given) 

7. P3 = 4.432 
U3 = 0.280 
83 = 0 . 930 

8. From = 1.288, Ms = I.966 and Ms = O.583 

ai 1 ^2 

as found in shock tables 
From equation (B1) U = -0.340 

Shock Wave Crossing Transverse Slits 


A schematic diagram of this interaction is given in sketch (b 4) 




Shock 

wave 


Transverse sUts 
/ /Entropy discontinuity 
/ 

Reflected 
disturbance 



X 

Sketch (b 4) 


As shown here, the resultant shock moves upstream and crosses the slits. 
The pressure difference between the channel and the surge chamber causes 
fluid to flow out of the channel through the transverse slits. This mass 
removal results in an immediate expansion which interacts with the shock 
wave and reduces its strength and velocity. The cnange in strength of 
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the shock wave introduces an entropy discontinuity into the fluid, as 
well as a small reflected disturbance. The method of solution for this 
interaction is given in the following outline form: 

Given : 


1 . All conditions (P, Q, u, a, t|) in sectors 1 and 5 * 

2 . Velocity is constant across entropy discontinuity: u^ = U3. 

3. Q and ^ are constant across the reflected disturbance: 

\ = and Q4 = Q5. It should be pointed out here that if this 
reflected disturbance should be a compression wave, then the 
invariance of and Q from sector 4 to 5 will not be strictly 
valid. Inasmuch as the disturbance is extremely weak, however, 
the variation of and Q is so small that, to the accuracy of 
this analysis, it cannot be detected. 

h. Dimensions of the transverse slits. 

Find: 


1 . All conditions (P, Q, u, a, t^) in sectors 2 , 3^ 3md 4. 

2 . The velocity of the shock after it has crossed the transverse slits. 
Auxiliary charts necessary (See Appendix A) : 


1. 

(P 2 -Pi)/ai 

vs (Q 2 -Qi)/ai 

for normal shock. 

2 . 


vs a^/ajL for 

normal shock. 

3. 

Steady- flow 

normal-shock tables . 


Procedure : 


1 . Assume (Q2 “Qi)/Q'1 • 

2 . Obtain (Ps-Pi)/^*! from chart. 

3. Compute P^^ ^2^ ^2* 

4 . Obtain ^2 = ^3 from chart. 

5. With the aid of equations (15) arid ( 20 ), compute P^, u^, and 

ag. Since the average values of the speed of sound and flow 
velocity must be used, an iterative procedure is necessary to 
perform this step. This iteration can be performed best by first 
using u = U2 and a = a^. From the values of and a^ thus 
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computed, new values of u and a can be computed. Repeat this 
process until sufficient accuracy is obtained. 

6. Compute (Q3 - Q4) from equation (AI8) replacing subscript 1 by 

subscript and subscript 2 by subscript 3* Use the given value 
of in computing the factor (Q4 + Pg) . 

7. Compare Q4( computed) with Qg( given) and repeat steps 1 through 6 

until these quantities are equal. 

8. Compute P4 from equation (AI7) replacing the subscripts as indi- 

cated in step 6. Use from step 7 to compute Ug, a^. 

9. Compute the speed of the shock after crossing the transverse slits 

using equation (B1) . 


A numerical example of this procedure is presented with l/8-inch 
transverse slits on both the top and bottom of the channel and the follow- 
ing conditions given in sectors 1 and 5: 


Pi = 5.384 

ui = 1.274 

% 

Qi = 2.836 

ai = 0.822 

7R 

P 5 = 5.571 

Ug = 0.280 

^5 

Qg = 5.011 

ag = 1.058 

7R 


= 0.212 


Following the specified procedirre, 

1. Assume (Q2-Qi)/ai = 2.368. 

Obtain (P2-Pi)/ax = 0.177 from chart. 

Compute P2 = 5.529 U2 = 0.372 

Qa = 4.734 ag = 1.031 

From a^/ai = 1.255 ( = nJyB. = O.165 


2 . 

3. 

4. 


By iterative procedure 
Qg - Q2 = 0.179 

Pg - P^ = -0.096 


5. 
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and 

P3 = 5A33 U3 = 0.235 

Qg = U.963 Eg = 1.040 

6. (^4 ^3)M = 0.01+7 (n4-'^3)Acp = 0.0047 

tanh ('''l 4 -'n 3 )/ 4 cp = 0.0047 

Q4 + P3 = 5-011 + 5-433 = 10.444 Use (given) 

(Q4 - Qg) = (P4 - Pg) = 0.049 

7- Q4 = 5-012 (computed) as compared to 
Q4 = 5.011 (given) 

8. P4 = 5.484 

= 0.235 
E4 = 1. 050 

9. Using ^ = 1 . 255 ^ obtain = I.857 and Me, = 0.604 from 

ai ^2 

shock tables. From equation (B 1 ) U = - 0.251 


Expansion Wave Overtaking Shock Wave 


The leading-edge shock wave of the strong compression pulse is fol- 
lov^ed by an extended expansion zone. This expansion zone propagates 
upstream more rapidly than the shock and^ hence^ overtakes and weakens 
the shock. The process is actually a continuous one but^ in applying 
the method of characteristics ^ the expansion zone is considered to be 
built up stepwise from a series of expansion waves. 


NACA TN 269k 


33 


A diagram of this process is shown in sketch (B 5 ) . The overtaking 
expansion wave interacts with the shock wave and reduces its strength 



and velocity. The change in strength of the shock wave introduces an 
entropy discontinuity into the fluid, as well as a small reflected dis- 
turbance. The method of solution for this interaction is identical to 
that presented previously for the interaction of two Q shock waves. 
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APPEiroiX C 

DETERMINATION OF THE INTIAL VALUES OF THE 
CHARACTERISTIC QUANTITIES 


The velocity profile for a given pulse has been superimposed upon 
the steady-flow velocity profile at the time of impending interaction 
of the pulse shock with the steady-flow normal shock. Since the condi- 
tions in the region upstream of the disturbance pulse are unchanged by 
the presence of the pulse, the speed of sound and all other properties 
of the fluid in this region are known from the steady-flow conditions. 
It remains necessary, therefore, to evaluate only the speed of sound 
throughout the region downstream of the pulse shock. With the velocity 
profile of the channel already known, the characteristic quantities can 
then be calculated from equation (9)* 


The following procedure can be used to determine the speed of sound 
immediately downstream of the pulse shock: Let sketch (Cl) represent 

the flow conditions in the neighborhood of the upstream edge of the pulse. 

IIIIIIII n ! I ! I n ! II ! I ! I ! ! n Superimpose a velocity -U on the 

flow field, thus converting to a 
coordinate system which represents 
velocities relative to the upstream 
edge of the pulse as shown in sketch 
(C2). The local speeds of sound, 
ag, and at>, as well as the ambient 
pressure, temperature, and entropy, 
are Independent of coordinate system 
and remain unchanged. With reference 
to sketch (C2), the following rela- 
tions can be obtained: 


u - — 


Ua ^ 



Sketch (Cl) 




Ug-U 

T) })) 1 1 nn t i / i 

1)1 II))} t 1 t 


Sketch (C2) 


Equations (Cl) can be combined to give 


Un - U 


= M, 


Ufa - U 

% 


= M. 


(Cl) 


^a - % 


= Msa - Msb $ 


(C2) 
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which expresses conditions across a stationary normal shock. Since ua^ 
Ujj, ag^, and pg are already known, the (quantities a^. Mg , and Mg can 

be determined from equation (C2) by a trial and error calculation with 
the use of any of the standard tables for flow across a normal shock. 

The solution to equation (C2) will also determine the pressure ratio 
across the shock from which the value of pq, can be computed. From the 
integrated form of equation (4), namely. 


An 


b 

a 



(C3) 


the value of *0^ can also be determined. 

To compute the speed of sound at all points dovnstream of the pulse 
shock, the following assumptions were made: 

1. The entropy is constant from the rear of the pulse shock, 
station b, to the end of the channel. 

2. The pressure at the end of the pulse, p^, is the same as 
the steady -flow pressure at that point. 

3 . The speed of sound along the pulse varies linearly with the 
flow velocity. 

From assumption (l), the expression for the entropy change between 
stations b and c is 


ATI 




(C4) 


By assumption (2) and equation (c 4 ), the speed of sound at the rear of 
the pulse, a^, can be determined since and are known. With the 
velocity profile of the channel already known, assumption (3) permits 
an evaluation of the speed of sound at all points along the pulse . The 
values of the characteristic quantities, P and Q, (eq. (9)) can now be 
evaluated at all points in the channel at the time of impending inter- 
section of the pulse shock with the steady-flow normal shock. 
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Figure 2. -Velocity profile of channel showing superimposed pulses in position of impending 
intersection with normal shack in stable equilibrium position. 
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Figure 3r interaction of critical and double critical pulse with steady 
flow normal shock for various amounts of mass removal. 
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(b) Velocity profile with length dimension corrected by scale 
factor 0.585 to produce desired pulse area. 
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Length along pulse tube, inches P 

(a) Characteristic diagram 


Figure 4.- Analysis to produce pulse with leading-edge velocity 
discontinuity, A~^, of -0.610 and velocity-profile area, 

(-^ vs inches) of 1.905. 
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Figure 5.~ Relation between puise tube chamber pressure and 

velocity-profile area. 
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Figure 6 .- Experimental apparatus 







Figure 7.- Pulse tube, 
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Figure 8- Design of puise tube, (aii dimensions in inches) 
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(b) Exposure, 1 microsecond. 


Figure 9.- Normal shock in channel under the influence of random disturbances 
originating downstream; transverse slits closed. . 
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(a) Exposure^ 1/25 second. 



(b) Exposure^ 1 microsecond « 


Figure 10 .- Normal shock in channel under the 
originating downstream; transverse 


influence of random disturbances 
slits open l /4 inch. 
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